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ABSTRACT 

We present ground-based mid-infrared imaging for 27 M-, S- and C-type 
Asymptotic Giant Branch (AGB) stars. The data are compared with those of 
the database available thanks to the IRAS, ISO, MSX and 2MASS catalogues. 
Our goal is to establish relations between the IR colors, the effective temperature 
T e ff, the luminosity L and the mass loss rate M, for improving the effectiveness 
of AGB modelling. Bolometric (absolute) magnitudes are obtained through dis- 
tance compilations, and by applying previously- derived bolometric corrections; 
the variability is also studied, using data accumulated since the IRAS epoch. 
The main results are: i) Values of L and M for C stars fit relations previously es- 
tablished by us, with Miras being on average more evolved and mass losing than 
Semiregulars. ii) Moderate IR excesses (as compared to evolutionary tracks) are 
found for S and M stars in our sample: they are confirmed to originate from 
the dusty circumstellar environment, iii) A larger reddening characterizes C-rich 
Miras and post-AGBs. In this case, part of the excess is due to AGB models over- 
estimating T e ff for C-stars, as a consequence of the lack of suitable molecular 
opacities. This has a large effect on the colors of C-rich sources and sometimes 
disentangling the photospheric and circumstellar contributions is difficult; better 
model atmospheres should be used in stellar evolutionary codes for C stars, iv) 
The presence of a long-term variability at mid-IR wavelengths seems to be lim- 
ited to sources with maximum emission in the 8-20 jum region, usually Mira 
variables (1/3 of our sample). Most Semiregular and post- AGB stars studied 
here remained remarkably constant in mid-IR over the last twenty years. 

Subject headings: Stars: mass-loss - Stars: AGB and post-AGB - Stars: carbon 
- Stars: Mira - Infrared: stars 
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Introduction 



Stars of low and intermediate mass (all those below M = 7-8 M Q ) terminate their 



Busso et al. 1999 ; 



evolution th rough the so-called Asymptotic Giant Branch (AGB) phase 
Busso 20071 ]. in which t hey lose mass efficiently thanks to stellar winds powered by radiation 



pressure on dust grains 



Habing 19961 ] . After this stage, they generate planeta ry nebulae and 



start a blue- ward path, which ultimately gives birth to a white dwarf [see also lHerwig 2005 ]. 



While moving along this path, AGB stars replenish the Interstellar Me dium with 



about 70% of all the matter returned after stellar evolution 



Sedlmavr 1994 1; th is is done 



through the formation of extended circumstellar envelopes [Winters et al. 20021 ] . Cool or 
cold dust, with its radiative emission in the infrared (IR), normally dominates the energy 
distribution of these sources, so that until recently the bolometric (apparent) magnitude 
of the most evolved AGB stars was difficult to derive, due to insufficient photometric 
coverage of the mid-far IR range of the electromagnetic spectrum. Traditionally, 



circumste 



Jura 1986 



l ar envelopes were studie d using fluxes provided by the IRAS satellite [see 



Jura fc Kleinmann 19891 ] . which had however a poor spatial resolution (hence 



some risk of contamination, especially at the longest wavelengths). Similarly, absolute 
magnitudes were very uncertain due to the difficulties in measuring the distance of these 
single, strongly variable stars. 

Only recently, the availability of a large IR database from space-borne telescopes like 
ISO, IRTS, MSX and the increased amount of ground-based mid-IR observations have 
substantially modified the situation. The quantitative studies of luminosities, colors and 
mass loss in the last evol utionary stages of mod e rately massive sta rs have consequently grown 



in number and quality 



Whitelock et al. 2006 



Feast et al. 2006|. In parallel, 



distances for AGB stars have been in part upgraded 



Bergeat fc Chevallier 20051 ] . The results 



Hipparcos 



may still include inconsistencies (Feast, private communication) but a full revision of the 
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Hipparcos catalogue is still waited for van Leeuwen 20051 ; Ivan Leeuwen &: Fantino 2005 ] . 



Recent works have therefore 
distances 



Menzies et al. 2006 : 



argefy explo ited Period-Luminosity relations to infer the 



surveys are becoming available 



east 20 07|] . Contem porarily, the first results of Spitzer's 



Ziilstra et al. 2006 ]. 



In the above modified scenario we started a reanalysis of galactic AGB stars at 
IR wavelengths, aimed at improving the determination of their energy distribution, 
mass loss and absolute magnitudes, making use of recent IR photometric data and of 
reliable distance estimates. With such works one can now finally compare homogeneous 
samples of mass losing stars in the Galaxy with similar data sets in Magellanic Clouds 



van Loon et al. 2001 



van Loon et al. 20051 ] . in order to study the dependence on 



metallicity of general stellar properties. 



In a previous paper 



Guandalini et al. 20061 . hereafter Paper I] we considered C stars 



as observed by the ISO-SWS and by the MSX experiments. Here we extend our work by 
analyzing a sample of 27 AGB stars of classes M, S and C (listed in Table 1) for which 
we obtained ground-based infrared imaging in the 10 /xm window. Most of the stars 
we observed were also the object of measurements by the above quoted space-borne IR 
telescopes, so that we can now compare and integrate results from different experiments 
and different epochs, as a check of the quality of the available infrared database and of the 
source variability. 

This paper is organized as follows: in section 2 we briefly discuss our IR camera, used 
for making the ground based observations. In section 3 we present the photometric data 
obtained through it at wavelengths longer than 8 //m, integrated by the near-IR archive 



observations of 2MASS [Cutri et al. 2003| and by the available estimates for mass loss 
rates, distances and a few other relevant parameters. In section 4 we then use the database 
thus constructed to derive bolometric magnitudes (both apparent and absolute), colors, and 
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correlations of photometric properties with mass loss rates. Then section 5 addresses the 
long-term variability issue, by going back to the IRAS catalogues in order to compare the 
available data over a time interval of about 20 years. Finally, in section 6 we outline some 
general conclusions, underlining the main remaining problems that remain to be solved for 
allowing a satisfactory match between photometric observations and stellar modelling. 



2. The Mid-IR Camera TIRCAM2 and its Calibration 



Our sample stars were observed between 2001 and 2004 at the Italian Infrared telescope 
(TIRGO), located at 3200 meters over the sea level, on top of the Gornergrat, in Switzerland. 
The data were collected using our mid infrared camera TIRC AM2 (Tirgo Infra red Camera 



version 2), an upgrade of a previously available instrument 



Persi et al. 1994 1. 



TIRCAM2 uses a Rockwell (then Boeing) high-flux Si:As BIB 128 x 128 array 
(HF-21) and is equipped with five narrow-band filters (10% bandwidth) between 8 and 
13 /im, with the N broad-band filter, and with a circular variable filter (CVF) having a 
spectral resolution of 3% in the 8-14 /xm range. The optics of the camera produce a plate 
scale at TIRGO of 0.77"/pix. The array, the optical system and the filters are assembled 
in a liquid-He-cooled dewar (HD-3[8]) of the Infrared Labs. Inc. (Tucson, Arizona). 
The readout electronics of the array, fully developed by our home Institutes, and the 
other general cha r acteristics of the inst rument have been previously described elsewhere 



Persi et al. 2002 



Corcione et al. 20031 ] . The absolute fluxes of some infrared standard 



stars in our photometric system are shown in Table 2. They have been derived from the 



spectral energy distributions (SEDs) published by Cohen and coworkers Cohen et al. 1999 



and references th erein], after a convolution with the spectral response of our filters 



Persi et al. 2002 1. 
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The narrow-band images at 8.8, 9.8, 11.7 and 12.5 /im of our sample of AGB stars 
were taken in the standard chopping & nodding technique, in order to remove the sky and 
telescope emission background. Images of standard stars were obtained during the nights 
at similar air mass as our targets, for flux calibration. They were also used to estimate 
the point-spread function (PSF) at each wavelength. The mean PSF (FWHM) obtained 
during our observing runs was ~ 3.2" in diameter at 8.8 /xm and ~ 3.8" at 12.5 /im. All the 
observed AGB stars appear as point-like sources at this spatial resolution. The on-source 
integration times were from 720sec for the faintest sources, down to 120sec for the brightest 
ones, in all the filters used. The derived detection limit is approximately 0.7 Jy (1 a) in 
300sec of integration time. The photometric measurements were extracted from the images 
of AGB stars through the DAOPHOT tool, within the IRAF data reduction package 



Stetson 19871 ]. using an aperture of 6". 



3. Infrared Colors and Mass Loss of AGB Sources 

The flux densities of the sources observed, as measured with the technique mentioned 
in the previous section, are given in Table 3, together with their la statistical errors. We 
collected, as a comparison, near infrared fluxes for our sources, as published in the 2MASS 
catalogue, together with observations by the MSX satellite (when available) for wavelengths 
longer than those of our filters. These data are presented in Table 4. Finally, distances, 
mass loss rates and wind velocities, as deduced by the literature, are presented in Table 5, 
together with the relevant references. 

By calibrating the data of Table 3 with the zero-magnitude fluxes deduced by standard 
stars in Table 2, we obtain the infrared colors of the sources. Examples of color-color 
diagrams are shown in Figure [1] and in Figure [2j Here TIRCAM2 sources are represented 
by big symbols: open-starred ones refer to M-type (oxygen-rich) stars, plain open ones to 
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Leech et al. 2004 1. 



S-type stars. Filled symbols indicate C-rich sources. Small dots refer to the sample of C 
stars we discussed in Paper I, with colors deduced from ISO and MSX observations, as a 
comparison (see Figure 6 in Paper I). We recall that the ISO data were taken from the las t 
on-line database available, where SWS measurements are fully calibrated 
From them we derived photometric estimates in the TIRCAM2 photometric system by 
convolving SWS spectra with the response curves of our filters (see Paper I). 

The bulk of our data is generally distributed along the line of black body emission in 
Figure [IJ while a larger dispersion is present in Figure [2J 

The spread in Figure [2] derives from the fact that, for sources that are presently on 
the AGB (Semiregulars or Miras), color indexes computed in the 8 — 14/im region are 
affected by the crowding of variable emission (or absorption) features characteristic of these 
wavelengths (e.g. at 11.2 — 11.7/im for C-stars, from PAH and SiC; and at 9.8 /im for 
O-rich sources, from silicates). When colors outside the 8 — 14/im band are used (Figured]) 
excess emission at 21 /im is seen for a number of post- AGB, C-rich TIRCAM2 sources (as 
also shown by ISO and MSX data, see Paper I). Emission at this wavelength was known 
;o be a distinctive property of a few post- AGB C stars and C-rich pre-planetary nebulae 



Kwok et al. 20021 ; IZhang et al. 20061 ] . However, so far the number of known sources with 



this feature was rather limited: in our data (both here and in Paper I) the 21 /im emission 
appears as a general property of post-AGB sources. The above simple considerations 
confirm the fact that mid-IR colors are effective tools for classifying AGB sources. When 
one deals with infrared stars, mid-IR colors identify C-rich AGBs through their color excess 
at 11.7/im, and C-rich post-AGBs through their 21 /im excess. 

Figure [3] then shows the relations linking the IR colors to mass loss. Again, the best 
measurements selected for C stars in Paper I are included for comparison. It is clear that 
our new C-star data follow the trends already established, as expected. It is however also 
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evident that the (few) O-rich sources of the sample behave differently, displaying lower 
mass-loss rates for a given value of the color, especially when this last makes use of fluxes in 
the 10/zm region (panel "a" in Figure [3]). Should this hint be confirmed, this would mean 
that the relation between dust opacities and mass loss rates for M and S stars has to be 
different than for C-stars. We are now verifying this in a dedicated work on a much larger 
database (Guandalini & Busso in preparation). 



4. Color-magnitude Diagrams and Effective Temperatures 



By applying to our database the bolometric corrections discussed in Paper I, we can 
easily derive apparent bolometric magnitudes of the C-ric h sources studied. For O-rich S 
and M stars we applied instead the bolometric correction 



Guandalini et al. 20071 ] : 



M bol - Af [8.8] = ax 3 + bx 2 + cx + d 



where x is the K-[8.8] color, and the coefficients are a=— 0.0211, b=0.0812, c=1.0658, 
d=2.3026 (with a correlation r 2 = 0.989). Correcting for the distance (see Table 5) yields 
the absolute bolometric magnitude, through which the absolute color-magnitude diagram 
can be plotted. Examples are presented in Figures 0] and [5j 

In Figure @] the dashed area r epresents the zone oc cupied by AGB photospheres 



according to current stellar models 



Straniero et al. 19971 ] . As the Figure shows, essentially 



all the stars of our sample, much like those of Paper I, are reddened with respect to 
the models, by amounts that are widely disperse. Especially for the reddest stars, the 
displacements are dominated by the known presence of dust in the circumstellar envelope. 
In this respect, an additional source of concern is however introduced by the poor knowledge 
of molecular opacities for C-rich atmospheres , and by the fact that they are usually 



neglected in stellar codes Marigo et al. 19981 ] . Due to this, it is known that the model 



- 9- 



effective temperatures of very cool, C-rich AGB models, in the region of thermal pulses, are 
unreliable and generally largely overestimated. For strongly variabl e, dynamically perturbed 



atmospheres the same physical meaning of T e ff becomes doubtful 



Uttenthaler et al. 2007 1- 



Uncertainties in model T e ff values have very different consequences on the colors of 
M-type and C-type stars. As an example, let's consider, for AGB stars of the two classes, 
the effects of a change by 0.1 dex in log T e ff. For M stars this corresponds e.g. to moving 
from type M3 to type M8: the (J-K) color difference in Figure H] is less than 0.2mag. On 
the contrary, using for C stars the classifications and the color calibrations by Bergeat et 



al. (2001), we see that the same s 
the class CV7 



lift, corresponding e.g. to moving from the class CV6 to 
Bergeat et al. 20021 ]. may imply a color difference (J-K) ~ 1 [see Figure 4 in 
Bergeat et al. 20011 ] . Similar considerations would hold for the [K-12.5] color in Figure 
We conclude that, while the colors of O-rich AGB stars are not significantly affected by 
uncertainties in the atmospheric models, for C-stars the situation is worse. At least the less 
reddened stars in Figure H] might show displacements from the model areas partly due to 
the real presence of dust, but partly also induced by errors in model T e ff values, hence in 
atmospheric opacities. 

A special word of caution must be added for the positions of post-AGB objects in the 
two plots discussed here. Sometimes they have a complex SEDs, related to the fact that 
the central star begins to be detached from the circumstellar shell and shines at a relatively 
high temperature, as typical of a yellow supergiant. In such cases the integral of the IR 
flux (and in general the procedure described in Paper I) is not sufficient to determine the 
bolometric magni tude properly, because the optical contribution is not negligible [see e.g. 



KwoketaL 1999 



Fig. 1]. For post-AGB objects we have therefore only lower limits to the 
bolometric magnitudes, which approximate the real values with an error that is different for 
different cases. 
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In order to illustrate the cautious remarks made above for the relations between the 
color excess and the presence of dust, we selected in our sample of TIRCAM2 sources, and 
in the "best set" of sources in Paper I (those with astrometric distances and very detailed 



SEDs up to beyond 40 um) a group of stars 



Qlofsson et al. 1993a 



Olofsson et al. 1993b 



or which independe nt estimates of T e tf exist 



Bergeat et al. 20011 ] or can be inferred from 



published color calibrations [Bessell et al. 19981 ] . For these objects we plot the absolute HR 
diagram (M& ; versus T e ff) in Figure [6) In the Figure, the tracks on the left refer to AGB 



models for various masses and metallicities 



Straniero et al. 20031 ] . It is evident that, while 



the observed magnitudes fall in the range expected by these calculations (which include 
minimal or no overshooting), a large number of sources, and in particular the C-rich ones, 
have temperatures much lower than those predicted: these last, lacking a treatment for 
C-rich molecular opacities, can be in error by up to 30%. 

The effects of absorption features from molecules like CO, CN and C2 was shown 
to be extremely sensitive to the effective temperature and of lar ge consequence fo r the 
near infrared colors by an hydrostatic analysis of a few C stars [Loidl et al. 20011 ] . The 
influence of molecular opacities for the changing composition of AGB envelopes gradually 
enriched by the TDU process, and in particular for C-rich sources, was also addressed by 
Marigo (2006). She showed how a large decrease in T e ff can be induced by the inclusion 
of opacities from molecules made of CNO and hydrogen. In models of about 2 M Q a shift 
in Log T e ff of almost 0.1 dex can be obtained (see Figure 6 in that paper). Just as an 
example, by reducing of this amount the temperatures of the model tracks in Figure [6] (in 
order to simulate the effect of molecular opacities) we get the curves at the right side of 
the plot (those with crosses superimposed). It is evident that a proper atmospheric model 
might be in many cases sufficient to yield the values of T e fj derived from observations. We 
warn that the new tracks must be seen only as a rough example, plotted for illustration 
purposes. In fact, Marigo (2006) computed a synthetic AGB evolution (while the tracks 
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in Figure [6] were derived from complete AGB models), and her assumptions 



Straniero et al. 2003| . What one 



or dredge-up 



were rather different than those in the models we used 
would really need here is a set of complete AGB models, including a proper treatment of 
surface opacities for the changing envelope composition during the TDU process, something 
that unfortunately does not exist yet. 

The above discussion says that, for C-stars, disentangling the atmospheric opacity 
effects from the color excess due to dust will become quantitatively meaningful only with 
large samples of good mid-IR observations, determining the properties of dust for sources 
where T e ff has been independently determined (e.g. from spectra). This in its turn would 
help stellar modelers to construct more detailed opacity tables, calibrated on observations. 



This is a relevant target to be pursued. In fact, the compilation o 
catalogues was started by several group s years ago [see for example 



Cioni et al. 2001 



years 



accurate infrared 



van Loon et al. 1999 



Bcrgeat et al. 2002 



Le Bertre et al. 2001 1 . It has c ontinued to be an essential tool in recent 



Le Bertre et al. 2003 



now [Bergeat fc Chevallier 200,4 iLe Bertre et al. 2005 



Cioni et al. 2003 1 and remain s important 



Whitelock et al. 2006I ]. It will 



also be one of the key projects of the Antarctic telescope IRAIT (International Robotic 
Antarctic Infrared Telescope) that we recently developed and that will be operative at 
the Italo- French base of Dome C s tarting from the 2007-2008 Antarctic campaign [see e.g. 



Tosti et al. 2006 



Busso et al. 2006 1 



5. Infrared Variability 

Figures [7J and M show the available information on the IR SEDs for two groups 
of sources in the list of Table 1. The plots include data from the IRAS-PSC (Point 
Source Catalogue), IRAS-LRS (Low Resolution Spectra), ISO and MSX, together with 
our TIRCAM2 measurements. Figure [7J shows distributions that, despite their different 
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appearance, share the property of being non- variable over a time interval of almost 20 years. 
The figure contains rather heterogeneous measurements: photometric data (from the IRAS 
Point Source Catalogue, from MSX and from our ground-based observations) are compared 
with spectroscopic information from ISO-SWS and IRAS-LRS. For our purposes this is 
sufficient: as verified in Paper I, the proper convolution of ISO-SWS and IRAS-LRS spectra 
with the response of our filters, yielding a homogeneous photometric database, would not 
affect the flux levels by more than 5%, which is well inside the internal uncertainty of each 
set of data used. A constancy of the IR energy distribution is the more common behavior 
displayed by our sources, being shared by exactly 2/3 of the AGB stars in our sample (18 
out of 27). This property characterizes stars that are very different from one another. They 
include sources with minimal IR excess (usually Semiregulars) , in which the SED is peaked 
in near-IR; but they also include evolved (post-AGB) objects, in which the maximum 
emission is at very long wavelengths (from 20 to more than 40 /im), due the dominant 
effects of cold, distant dust. 

In contrast, Figure El shows the behavior of intermediate objects, usually Mira variables, 
in which the emission peaks near 10 /xm (and sometimes is rather flat up to about 20/xm), 
efficiently powering the typical features there present for O-rich and C-rich dust. Nine 
sources share this behavior (RAFGL sources n. 190, 809, 865 and 954; IRC+60144; RU 
Vir; CIT 6; CW Leo; IRC+40156). It seems, therefore, that long-term mid-IR variability is 
not a common property of AGB stars, but is restricted to a special class of sources, in the 
special evolutionary stage when most of the flux is re-radiated by circumstellar layers of 
one to a few hundred K. Here stellar pulsation (which is of large amplitude in the optical 
bands) should effectively transfer energy to dust shells, which are quite opaque down to the 
mentioned temperatures. One can guess that this is the phase where the coupling between 
the photosphere and the circumstellar envelope is most efficient: for bluer sources, not 
enough dust is created, probably by a mass loss rate which is still not strong enough. For 
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redder objects, circumstellar dust becomes cold and distant, probably detaching itself from 
the central star. 

One might a priori argue that the Semiregular sources might look stable only because 
of their small IR fluxes: when the flux is very low any variability might be more difficult 
to disentangle from the background noise. In order to make our suggestions more secure 
we looked for sources (both variable and non-variable) for which the ISO-SWS instrument 
offered repeated observations. We found a few interesting cases. Examples are presented 
in Figure [9] for the C-rich semiregular star S Scl and for the C-rich Mira variable V Cyg. 
Even looking, in the figure, only at homogeneous data sets (those from ISO, covering about 
2 years) our previous suggestions seem to be confirmed. It is indeed clear that the Mira 
variable does vary as we mentioned, and it can also be understood how the non-variability 
of S Scl in mid-infrared in the time interval covered by the ISO data (at a level better than 
10%) does not descend from an insufficient precision, since the minimum flux is of the order 
of tens of Jy. 

The idea that the variable sources are those for which radiation pressure on dust 
grains powers mass loss suggests itself, though we can look at this only as at a tentative 
interpretation, because the source statistics and also the number of available points per 
source are rather limited. This hypothesis deserves now to be verified through modelling 
of the radiation transfer in the dust-condensing envelope, also to determine whether IR 
variability corresponds to a specific evolutionary stage, or can be encountered repeatedly 
(the same dilemma presented by optical Mira variability). 
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6. Conclusions 

The data presented in this note confirm the usefulness of mid-IR colors taken from 
ground based telescopes, in describing the properties of mass-losing AGB stars. Our 
photometric study yields results compatible with those of Paper I, and begins to extend the 
analysis to a number of O-rich (S and M) sources. These will now be the object of dedicated 
papers, based on a large space-borne observational database similar to that of Paper I. 

In general we showed how, for most semiregular variables and post-AGB stars, the IR 
fluxes longward of a few microns are insensitive to the variations in the stellar photosphere, 
and remain essentially constant in time over rather long time intervals (tens of years). 
This however does not apply to Mira variables, and in general to AGB sources having 
their maximum emission near 10 /im. Such objects show remarkable variability in the 
emission/absorption features (with changes from emission to absorption and back) and/or 
in the global flux, showing that warm dust quickly reacts to the large-amplitude surface 
pulsations, with large changes in concentration and temperature. This might suggest that 
dust-driven winds are mainly associated to the Mira stage, a hypothesis to be further 
verified. 

Our results also point out how urgent it is to match stellar evolutionary codes with 
reliable model atmospheres, especially for carbon stars, including molecular opacities. This 
would allow one to predict reliably the effective temperature of C-rich atmospheres and 
hence to give a quantitative meaning to the IR color excess in terms of dust emission, 
without remaining uncertainties from poorly understood photospheric opacities. 
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FIGURE CAPTIONS 

Figura 1 An example of a color-color diagram for the observed stars. Big symbols are from 
TIRCAM2, small ones from Paper I. See text for comments. 

Figura 2 Another example of a color-color diagram for the observed stars. Again, big symbols 
are from TIRCAM2, small ones from Paper I. See text for explanations. 

Figura 3 Mass loss rates as a function of the infrared colors for the sample of stars in this 
paper (large dots) and for the best-measured sources of Paper I (small dots). As usual, open 
symbols refer to S stars, filled symbols to C stars. 

Figura 4 The color-magnitude diagram for our sources, using the J-K color as a temperature 
indicator (see text for explanations). Small symbols are from Paper I. The "C-star limit" 
line refers to 2 M models of solar metallicity. 

Figura 5 Same as above, but using the K-[12.5] color in abscissa. This color samples dust 
of different temperatures (and also, in K, the residual photospheric flux); the horizontal 
extension of the data is a measure of the infrared excess due to the dusty envelope. 

Figura 6 The HR diagram of the few TIRCAM2 C-rich (full large dots) or S-type (open 
large symbols) stars for which we could find in the literatu re independent estim ates of T e ff. 



For S-stars T e jf is an extrapolation [from the relations by lBessell et al. 19981 ] . Small sym- 
bols refer to those C-rich stars from Paper I having astrometric distances and values of 
T P ff. The curves at t he right are from the AGB models at various masses and metallicities 



Straniero et al. 20031 ] . It is evident that, on average, C-rich sources show large temperature 



discrepancies compared to models sequences. The curves on t he right side are shifted by 



O.ldex in Log T e ff, in order to mimic molecular opacity effects [Marigo 2007 1 . 



Figura 7 The Spectral Energy Distributions of a few non-variable sources, as available from 
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our TIRCAM2 observations and from the IRAS PSC, IRAS LRS, ISO-SWS, MSX, and 
2MASS catalogues. SEDs with maximum emission in the near-infrared, as well as those 
with maximum emission longward of 20/im all show a constant flux in mid-IR. 

Figura 8 Spectral Energy Distributions of sources that show significant variability over the 
time elapsed from the IRAS to the TIRCAM2 observations. Again, data obtained with 
TIRCAM2 and data available from the IRAS PSC, IRAS LRS, ISO-SWS, MSX, and 2MASS 
catalogues are included. Only sources for which the emission is maximum in the range 
8 — 20/im do appear to be variable, independently on their composition (C-rich or O-rich). 

Figura 9 ISO-SWS Spectral Energy Distributions of two Carbon Stars: a Mira variable (top 
panel) and a Semiregular (bottom panel); they were observed repeatedly by ISO during 
about 2 years. On this rather short time scale, but with homogeneous data, our suggestions 
about mid-IR variability of Miras and flux constancy of Semiregulars seem to be confirmed. 
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Tabella 1 The TIRCAM2 Sample of AGB Stars. 



Number 


IRAS Name 


Other Name 


Coordinates 


Var. Type 


Type 


1 


01144+6658 


RAFGL 190 


01 17 51.62 +67 13 55.4 


P 


C 


2 


03186+7016 


RAFGL 482 


03 23 36.57 +70 27 07.5 


M 


c 


3 


04307+6210 


IRC +60144 


04 35 17.45 +62 16 23.3 


S 


c 


4 


04395+3601 


RAFGL 618 


04 42 53.67 +36 06 53.2 


P 


c 


5 


04530+4427 


RAFGL 6319S 


04 56 43.28 +44 32 41.6 


- 


c 


6 


05405+3240 


RAFGL 809 


05 43 49.78 +32 42 06.8 


M 


c 


7 


05426+2040 


Y Tau 


05 45 39.41 +20 41 42.1 


S 


c 


8 


06012+0726 


RAFGL 865 


06 03 59.84 +07 25 54.4 


M 


c 


9 


06176-1036 


Red Rectangle 


06 19 58.22 -10 38 14.7 


P 


c 


10 


06291+4319 


RAFGL 954 


06 32 41.93 +43 17 15.3 


I 


c 


11 


07134+1005 


HD 56126 


07 16 10.26 +09 59 48.0 


P 


c 


12 


09452+1330 


CW Leo 


09 47 57.38 +13 16 43.7 


M 


c 


13 


10131+3049 


CIT 6 


10 16 02.27 +30 34 18.6 


S 


c 


14 


12427+4542 


Y CVn 


12 45 07.83 +45 26 24.9 


S 


c 


15 


12447+0425 


RU Vir 


12 47 18.41 +04 08 41.4 


M 


c 


16 


12544+6615 


RY Dra 


12 56 25.91 +65 59 39.8 


S 


c 


17 


06331+1415 


DY Gem 


06 35 57.81 +14 12 46.1 


S 


s 


18 


09076+3110 


RS Cnc 


09 10 38.80 +30 57 47.3 


s 


s 


19 


12417+6121 


S UMa 


12 43 56.68 +61 05 35.5 


M 


s 


20 


15492+4837 


ST Her 


15 50 46.62 +48 28 58.9 


s 


s 


21 


06297+4045 


IRC +40156 


06 33 15.75 +40 42 50.9 


p 


M 


22 


12277+0441 


BK Vir 


12 30 21.01 +04 24 59.2 


s 


M 


23 


13001+0527 


RT Vir 


13 02 37.98 +05 11 08.4 


s 


M 


24 


14059+4405 


BY Boo 


14 07 55.76 +43 51 16.0 


I 


M 


25 


14219+2555 


RX Boo 


14 24 11.63 +25 42 13.4 


s 


M 


26 


14371+3245 


RV Boo 


14 39 15.86 +32 32 22.3 


s 


M 


27 


16269+4159 


g Her 


16 28 38.55 +41 52 54.0 


s 


M 
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Tabella 2 Flux Densities of standard stars. 



Standard 


[8.8] 


[9.8] 


[11.7] 


[12.5] 




(Jy) 


(Jy) 


(Jy) 


(Jy) 


a Lyr 


49.69 


40.44 


28.48 


25.05 


f5 Gem 


152.74 


120.90 


88.67 


76.77 


a Boo 


883.72 


745.32 


524.94 


459.01 


(3 And 


306.36 


263.52 


200.79 


174.55 


/3Peg 


431.12 


376.00 


279.25 


249.01 


a Tau 


752.09 


646.86 


481.44 


419.57 
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Tabella 3 Our TIRCAM2 observations in the 10/im window. 



Source Number 


Epoch 


F[8.8] (Jy) 


F[9.8] (Jy) 


F[11.7] (Jy) 


F[12.5] (Jy) 


1 


16.1.2003 


73.5(5.4) 


29.2(8.2) 


78.9(11.3) 


109.6(14.4) 


2 


16.1.2003 


147.1(7.9) 


- 


146.8(17.6) 


154.6(19.1) 


3 


16.1.2003 


174.8(8.8) 


156.5(10.9) 


133.1(16.4) 


145.4(19.2) 


4 


6.12.2003 


254.6(27.1) 


380.2(42.8) 


490.2(26.7) 


488.0(28.4) 


5 


13.2.2002 


87.2(19.1) 


83.0(8.3) 


100.1(10.0) 


87.8(9.0) 


6 


16.1.2003 


185.8(9.3) 


158.2(11.9) 


143.1(17.5) 


169.5(21.5) 


7 


12.2.2002 


122.6(13.5) 


116.8(7.4) 


133.0(6.7) 


80.9(3.6) 


8 


6.12.2003 


364.9(28.9) 


414.1(41.1) 


477.9(34.6) 


379.1(43.2) 


9 


7.12.2003 


301.5(34.2) 


307.5(50.4) 


387.2(35.9) 


387.4(49.0) 


10 


15.1.2003 


63.4(6.6) 


66.2(14.6) 


129.1(12.6) 


81.6(24.) 


11 


6.2.2004 


5.6(1.8) 


10.7(2.4) 


7.3(2.8) 


20.6(2.7) 


12 


14.1.2001 


30255.(605.) 


- 


41782.(342.) 


38978.(667.9) 


13 


15.1.2003 


1789.(156.) 


1976.(190.) 


2226.(400.) 




14 


6.12.2003 


342.8(20.3) 


337.0(29.0) 


108.9(22.0) 


220.7(24.6) 


15 


16.1.2003 


253.0(27.6) 


252.4(24.3) 


245.7(39.3) 


208.8(25.5) 


16 


16.1.2003 


140.1(15.7) 


83.0(10.7) 


84.4(13.8) 


104.4(13.4) 


17 


6.2.2004 


5.3(1.6) 


16.4(3.3) 


34.7(10.5) 


24.5(4.5) 


18 


14.1.2001 


398.4(40.0) 


550.0(82.5) 


453.1(63.4) 


378.3(57.0) 


5? 


11.2.2002 


568.5(47.8) 


790.7(29.8) 


497.2(13.0) 


452.1(33.0) 


5? 


6.12.2003 


570.0(30.0) 


739.1(49.0) 


530.0(37.0) 


477.1(28.0) 


19 


11.2.2002 


4.5(0.5) 


- 


2.8(0.5) 


3.1(0.7) 


20 


4.2.2004 


155.7(8.0) 


205.2(33.0) 


201.4(50.0) 


215.5(28.0) 


21 


15.1.2003 


209.4(10.8) 


334.1(37.8) 


321.4(69.0) 


134.9(10.0) 


22 


4.2.2004 


224.6(11.3) 


217.6(32.6) 


118.4(29.6) 


211.6(27.5) 


23 


14.1.2001 


365.7(16.4) 


496.4(24.0) 


394.3(22.4) 


432.5(54.3) 


24 


12.2.2002 


77.7(8.6) 


72.2(1.3) 


51.6(2.6) 


40.2(2.3) 


25 


11.2.2002 


694.8(58.4) 


922.2(34.7) 


795.8(20.5) 


729.2(52.9) 


26 


6.2.2004 


106.7(13.9) 


151.3(47.2) 


75.1(21.0) 


86.4(8.6) 


27 


12.2.2002 


461.6(50.9) 


427.4(27.0) 


340.4(17.3) 


293.6(13.2) 
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Tabella 4 Source fluxes in the 2MASS, MSX fc ISO filters (Jy) 



Source no. 


Var. type 


J 


H 


K 


14,6 


21,3 


Data orig. 


1 


P 


0.01 


0.02 


0.01 


96.3 


110 


ISO 


2 


M 


0.09 


1.1 


6.8 








3 


S 


4.3 


16.0 


48.9 








4 


P 


< 0.01 


0.03 


0.20 


651 


1260 


ISO 


5 




< 0.01 


0.03 


0.35 


62.0 


51.2 


MSX 


6 


M 


0.05 


0.53 


3.6 


227.5 


147.8 


MSX 


7 


S 


242 


495 


481 


59.1 


36.4 


MSX 


8 


M 


< 0.01 


0.20 


2.5 








9 


P 


3.7 


9.0 


23.0 


377 


472 


ISO 


10 


I 


8.1 


22.7 


39.3 








11 


P 


2.9 


2.1 


1.5 


42.0 


116 


ISO 


12 


M 


2.7 


74.9 


469 


26309 


18194 


ISO 


13 


S 


2.5 


17.1 


91.4 








14 


S 


641 


1331 


1316 


53.8 


33.5 


ISO 


15 


M 


9.7 


27.2 


82.2 


91.6 


54.3 


ISO 


16 


S 


218 


425 


464 


50.4 


24.7 


ISO 


17 


S 


85.4 


156 


143 


12.7 


7.8 


MSX 


18 


s 


3065 


4324 


3743 








19 


M 


26.3 


43.4 


41.4 








20 


s 


804 


1162 


1098 


149 


104 


ISO 


21 


p 


11.5 


33.8 


78.3 








22 


s 


966 


1448 


1306 








23 


s 


1145 


1843 


1770 




167 


ISO 


24 


I 


630 


1033 


826 








25 


s 




4210 


3948 


584 


422 


ISO 


26 


s 


512 


785 


700 


76 


66 


ISO 


27 


s 


3841 


5627 


4760 


273 


150 


ISO 
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Tabella 5 Relevant parameters and distances of the sample stars. 



Source no. 


Var. type 


M [M Q /yr] 


Type (ref.) 


v [km/sj 


in i 

d |kpc| 


Type ref. 


1 


P 


2.45E-05 


L 


18.2 


2.79 


G 


2 


M 


1.05E-05 


L 


12.2 


1.97 


G 


3 


S 


6.33E-06 


L 


18.5 


1.03 


G 


4 


P 


2.00E-04 


Meil998 


19.5 


1.70 


Meil998 


5 




1.46E-05 


L 


20.2 


2.60 


G 


6 


M 


2.40E-05 


L 


28.0 


2.01 


G 


7 


S 


1.60E-06 


B 


11.0 


0.74 


B 


8 


M 


1.16E-05 


L 


15.8 


1.47 


G 


9 


P 


1.00E-04 


Men2002 


5.0 


0.71 


Men2002 


10 


I 


6.37E-06 


L 


21.4 


2.19 


G 


11 


P 


2.23E-05 


L 


10.7 


2.40 


Hon2003 


12 


M 


3.30E-05 


B 


14.7 


0.15 


B 


13 


S 


6.50E-06 


B 


17.0 


0.41 


B 


14 


S 


1.40E-07 


B 


8.5 


0.26 


B 


15 


M 


2.30E-06 


B 


18.4 


0.68 


B 


16 


O 

o 


4.40E-07 


B 


1 A A 

10.0 


0.55 


B 


17 


S 


6.27E-08 


G-dJ 


8.0 


0.56 


G-dJ 


18 


s 


3.92— 5.20E-08 


G-dJ 


7.2 


0.12 


Hip 


19 


M 








1.15 


G-dJ 


c\r\ 

20 


S 


4.52— 7.21E-08 


G-dJ 


9.1 


0.31 


TT* 

Hip 


21 


P 


1.00E-05 


L 


16.3 


1.60 


L 


22 


s 


5.49E-07 


W 


7.5 


0.18 


Hip 


23 


s 


3.70E-07 


L 


9.3 


0.14 


Hip 


24 


I 








0.14 


Hip 


25 


s 


6.48E-07 


L 


10.2 


0.16 


Hip 


26 


s 


5.49E-07 


L 


8.1 


0.39 


Hip 


27 


s 


1.43E-07 


L 


8.5 


0.11 


Hip 



